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On thermal phase structure of deformed Gross-Neveu model 
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Abstract 

We illustrate the phase structure of a deformed two-dimensional Gross-Neveu 
model which is defined by undeformed field contents plus deformed Pauli matrices. 
This deformation is based on two motives to find a more general polymer model 
and to estimate how g-deformed field theory affects on its effective potential. There 
found some regions where chiral symmetry breaking and restoration take place 
repeatedly as temperature increasing. 
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Gross-Neveu (GN) model and its relevant models have been studied intensively as 
a simple, but clear, model of chiral symmetry breaking in various situations; under the 
circumstances of thermal 0, electromagnetic H |3[] or curvature 0] backgrounds (and 
references are therein). 

Another plausible application of these (two-dimensional) models is an effective theory 
to one-dimensional systems such as a polymer (for a review see ||). It is shown in 
that the N = 2 GN model is equivalent to a continuum model of polyacetylene. Further, 
the first order phase transition of polyacetylene from solitonic (chiral broken) phase to 
metallic (symmetric) phase of is formulated as a response to the doping concentration 
(chemical potential), i.e. the number of dopants per carbon atom, using large N ap- 
proximation J7|. The critical value of chemical potential coincides with an experimental 
value (Recently, the 1/N correction to this model has been calculated ||). There is 
also a discretized version (fermionic vector model H) which has been studied as one of 
polymer models on discretized surfaces (random polymers) or as string theory and lower 
dimensional quantum gravity [ |T0|j . 

As shall be shown later, we immediately see that the effective potential of the two- 
dimensional GN model (defined in the large N leading order) under a static background 
corresponds to the electrostatic potential of a specific charge distribution on a half-line. 
Consequently, solitonic-metallic phase transitions of a polymer can be regarded as phase 
transitions concerning its electrostatic states as well. This suggests that there exist 
polymer models as many as electrostatic potentials. On the one hand, one may expect a 
universal structure of electric charge distributions on a polymer. On the other hand, we 
have only a few GN models coupled to background fields, which play the role of charge 
distributions. Hence, we introduce a new example of electrostatic potential, which is 
under an oscillatory distribution, deforming the effective potential of the GN model 
coupled to a constant background gauge field. This approach would broaden the range of 
application of model within the method of large N calculation. Since we deform the Pauli 
matrices, our deformation may correspond to a kind of fluctuating noise, quon particles 



II] or a violation effect similar to the Pauli principle violation |L2 , which are formulated 



using g-deformation. 



From a viewpoint of high energy physics, the deformation possesses the following 
meaning. Motivated by unusual spacetime structure like non- commutative geometry 



13], which would be anticipated in an extremely high energy scale, there are several 



g-deformed field theory constructions We would like to focus on a pure effect from 
g-deformation rather than from gravitational effect here. However, all of these deformed 
models are not yet viable to calculate effective actions performing path integrals be- 
cause of their mathematical intricacy — fermions belong to quantum Lorentz group, and 
deformed Pauli matrices are defined by multiplications of g-factors which are governed 
by an appropriate quantum R- matrix fl5fl . In order to catch a glimpse of g-deformed 
field theory aspect, we employ a likely approximation which definitely makes calculation 
possible. Namely, we will use ordinary (undeformed) fermions plus a common g-factor 
(g = e em ) th the Pauli matrices. 

In this report, we investigate the thermodynamics and its phase structure of the 
deformed GN model. If we take the vanishing limit of an electric background (or of the 
deformation parameter), our model is exactly reduced to the thermal generalization of 
the original model discussed in H. The phase diagrams we shall show here cover every 
known phase boundary of the undeformed GN model accordingly. One of our interests 
is thereby how the only 1st order transition point extends into a finite region of phase 
space. 

Our starting model Lagrangian (before deforming) is 

C = fa-fid^ - ieAjip - —a 2 - uijjip, (1) 

where N stands for the number of flavors, and summations over flavors are implicit. 
is a constant background field. The large N effective (bare, D = 2) potential of this 
model is expressed in the following form 

1 roc I 9 

V(a; - V(0, = Yy o 2 + J q ds—Q{s£,)e- s " , (2) 

where £ is a background constant as we shall see later. The second term (quantum 
correction) on r.h.s. of (Q) can be interpreted as the electrostatic potential of an effective 
charge (form factor) Q with damping factor e _so " 2 on the half-line parametrized by s, while 
the first (classical) term is a constant in a sense. Thus, the function Qe~ sa2 represents a 



charge distribution in one dimension, and a~ 2 corresponds to an effective range of electric 
force. In the metallic (chiral symmetric) phase, < a >= 0, the effective range becomes 
infinity, and the potential behaves as the Coulomb potential of charge Q. £ plays the role 
of einbein as it stands in (|2|), but we can see its meaning more clearly in other models. 
For example, in the higher dimensional (magnetic dominant) cases of ([!]), Q and £ are 
given by @ 



Q(s) = scoth(s), f = eV-E 2 + B 2 , (3) 

where E is assumed to be small and perpendicular to B. When E — 0, £ _1 is called 
the magnetic length which represents the radius of cyclotron motion of each electron 
in classical mechanical picture, and hence electrons are spaced by In the case of 
constant curvature background, we have, for example fT7| , 
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where a weak curvature approximation is assumed. In each case, £ _1 represents a specific 
length of the system and plays the role of einbein in the proper-time integral. Owing to 
this property, we may suppose a lattice system of spacing and this fact might have 
a relation to carbon lattice systems. 

Let us go back to the D = 2 story. Obviously, in chiral symmetry breaking the- 
ories, there is no sense of taking D = 2 limit in (0) because of no magnetic field in 
two dimensions. However, as mentioned previously, it still makes sense as a problem 
of 1-dimensional electrostatics and generalization of polymer model, since what func- 
tional form of Q we choose corresponds to a definition (deformation) of polymer model. 
Furthermore, from a field theoretical viewpoint, D = 2 is just a toy model which sim- 
plifies most of equations, however dimensionality is not a crucial factor for the purpose 
of observing a pure effect of g-deformation. One can repeat the same analysis in higher 
dimensions. We therefore define every quantity of D = 2 theories in the formal limit 
D -> 2. 

In order to study an oscillatory function, we further deform (||), and our choice of Q 

is 



Let us briefly see why the choice is a deformation of (|3|). In the process of deriving 
(H), using the dimensional regularization, we have introduced the following slight phase 
difference between the gamma matrices and the auxiliary field a 



o - a, lfl - e e7r V (6) 

This procedure brings a slight violation of chiral symmetry (and of course the replace- 
ments of 7 M in ([!])). This gives a naive g-deformation (q = e em ) of the gamma matrices 
in the D = 2 model. Since Q of eq.(^|) is originally proportional to 

Trexp^esF^] x exp[-Tr [in ], (7) 

~~ \ / m^ 

the phase deformation (|6]) is relevant to only the first trace (over the gamma matrices). 
Here we extract particular two cases concerning the first trace, 

cosh(es£) when e = n/2, (n E Z) 
cos(es£) when e = (2n + l)/4. 

Nothing happens to the charge distribution function in the former case, while the latter 
case gives rise to an oscillatory function without any singularity. When e is an integer, 
which corresponds to the undeformed model, we reproduce the theory (|3[). This case will 
be reported in a forthcoming paper. Choosing the latter CcLSG clS cL deformed GN model, 
we arrive at (|5]). We should consider £ as a specific electric length under the charge 
distribution given by (|5p. 

Now, we are in a position to discuss the thermodynamic potential of our model defined 
by (0) with eq.(|5|) which is deformed from the theory (|3|) and is defined in the formal limit 
D — > 2 as already mentioned. The thermodynamic potential is obtained in the following 
form applying the method of [17|| to our model (notation and derivation are completely 



parallel to frT| ) 



1 2 tr[l] /■°°rfs@2(s 
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where T = and we have introduced the dimensional regularization and tr[l] means 
the trace of gamma matrix unit. The polymer model of [0] is a particular case of this 
potential (£ = 0, (3 — * oo and of course D — 2), and the case of £ = with a finite (3 is the 



thermodynamic generalization of the polymer model. Note also that the replacements of 
Q with (D and (|) reprod uce the thermodynamic potentials discussed in fl6| and [|l7j 
respectively. Imposing the condition 

r^&^^^L^Tn- (10) 
the renormalization of coupling constant A# is given by 

and the renormalized effective potential is therefore 



x [ ke-°> - l)^ e ^e 2 («^, isp + <xV'(l - 2s) ]. (12) 

In principle, phase transitions are classified by behaviour of potential minima, i.e., by the 
gap equation dV/da = 0, 

(We hereafter adopt the following value of the renormalized coupling constant 

which means the broken phase with the dynamical mass a — 1 at T — jj, — £ — 0.) 

However, there are more convenient equations which can be derived from the gap 
equation. These are in order: 
(i) 2nd order critical surface satisfies 



= r ^^lQ( s O{-e^^e2(s 2 J, isp + e-(l - 2,)} + 2se~% (15) 
defined by 



limA%(,^) = 0. (16) 
(ii) 3rd order critical line is determined by simultaneous solution of ( p!5| ) and 



f°° ds sa^t ^\v4vrs^ , 2/i . 4tt 
°=' j4^ Se Q{s() —^ S J'' S lP h (17) 



which follows from 



(iii) 4~th critical point is defined by simultaneous solution of fll5|) , ( |i~7|) and 

Hm(^)V(<7;£,/3,/x) = 0. (19) 
This is alternatively written in the following, but it has no solution in our case: 

oo 



/ ds n j v^s^ . 2u . 47i\ 
°=/„ (4^ SeW(?(s «— e »('f'V ) - <20) 



Remember that every solution of fllBD is not always on the true critical surface. 3rd order 
critical line cuts off irrelevant solutions of fll5|) away from the critical surface. Instead, 
1st order critical surface should be found from simultaneous solution of V = and 
dV/da = 0. 

Note that eqs.(|l|) and (|17|) are reduced to the following simple equations in some 
regions: 

(a) Limit of £ — > (T-/x plane): 

P D - 2 r(l - j ) = 2 r(^)ReC(3 - A i + i^-) for (15), (21) 

and 

ReC(5-D,i + A = for (17), (22) 
where £ is the generalized zeta function. These equations are exactly the same as those 



derived in [|18|] where various critical equations are reproduced from these equations. 

(b) Kmit ofT — > (£-axis, /i = 0): 

arc^e 1 -* + iMr 1 ) + -MC 1 ) = o for (is), (23) 

^-Him Mzrl) =0 for (17), (24) 

where ipo is 

/•OO 

^(«) = / ^ S - D / 2 (e— - l)Q(s). (25) 

JO 

(c) Leading order in weak \i expansion: 

= 1 ^[Q( S (){e- s (l-2,)-e s '' 2 (l + 2Vr 1/2 } + 2^] for (15), (26) 



= J o j^yjf^ (1 + 2 Vr V2 Q«) for (17). (27) 

The results on our phase structure determined by these equations in D = 2 are 
summarized in Figs. 1-6. In Figs. 1-5, we show various sections (£-T planes) of the critical 
surface from \i = to 0.7. In Fig.6, we show the solutions of ( |26| ) and (|27]). S and 
M express solitonic/metallic (in the original words, broken/symmetric) phases. On the 
solid/dashed lines, 2nd/lst order phase transitions take place. The points A, Ax, A 2 , 
B, Bi, B 2 are the tri-critical points which, strictly speaking, form two lines as explained 
below. 

Let us go into details of how the critical surface looks like in order of values of \i. 
For the interval < \i < 0.59, all the phase transitions are second order transitions 
which satisfy flTBD . In Fig.l, we show a representative diagram (/i = 0) for this interval 
of fi. Its critical values of £ and T are £ c = 0.75 and T c = 0.57. Similar diagrams 
as Fig.l appear, as fi increasing whereas T c decreasing. When fi = 0.59 (Fig. 2), we 
observe the point A which is just a turning point of the tricritical line A\-A-A 2 (see 
Figs. 2-5). This point is not 4-th critical because (^) has no solution as mentioned 
previously. The coordinates (£,T) of A are (0.72,0.08). All the transitions of Fig. 2 
except for the point A are second order transitions. When fi = 0.601 (Fig.3), another 
tricritical point B = (0.23, 0.3) appears just similarly as the point A. The point A now 
splits into A\ and A 2 . The dashed line A\-A 2 of 1st order transitions are determined 
by solving the equation (|i~2"l)= with (|T4|). In Fig. 4 (fi = 0.607), we figure out two 
tricritical lines A\-A-A 2 and B\-B-B 2 (The dashed lines A\-A 2 and B\-B 2 are of course 1st 
order). B u B 2 , A u A 2 are (0.1, 0.31), (0.31, 0.285), (0.74, 0.094), (0.71, 0.07) respectively. 
B\ disappears into the negative region of £ at fi = 0.608 and T = 0.32. This is the 
well-known tricritical point [|TJ. Finally, in Fig. 5, the coordinates of B 2 , A\, A 2 of = 0.7 
are (0.55, 0.24), (0.8, 0.108), (0.69, 0.059) respectively. 

We show in Fig.6 a plane diagram as well, solving fl26|) and ([27]). In this case, we 
find a tricritical point at (/x, £) = (0.607,0.695) and 1st order transitions in fi > 0.607. 
This might mean that the approximation becomes invalid around this point (numerical 
analysis does not work well for some higher value of /i). This is a different feature from 
fj~7| which shows only second order transitions. 



In conclusion, let us summarize our study. We have found, taking a particular charge 
function, the 1st order transitions in certain regions of temperature and the electric length 
parameter £, whose small value region is continuously connected to the known transition 
point shown in 0. This smooth connection of 1st order transitions would seem to be 
natural. On the other hand, appearance of the tricritical line Ai-A-A 2 (and the concave 
structure of A1-A2) is a main difference between our model and the other models @ and 



17]. This difference comes out of oscillatory nature of the Q function. Although we 



have observed respective £-T sections of the critical surface, we could have also shown 
T-fi sections at several £ values. According to Fig. 6, assuming the first order line to be 
connected to the origin of Fig.5, the critical value of // decreases as £ increasing. This 
feature can be a replacement of other models. 

For example, the large N £ = model gives a perfect agreement with an experimental 
data and the 1/N correction reduces the critical chemical potential by 20%, which is still 
acceptable within experimental limits ||. Our results of critical chemical potentials for 
< £ < 0.695 stay within this 20% reduction and therefore experimentally acceptable. 
We can offer various models of polyacetylene in conformity with the choice of £. 

The concave structure of A\-A 2 curves shows that the phases around there sensitively 
respond to temperature, and phase transitions (mass generations) repeat several times 
as T increasing. This feature can not be seen in the undeformed cases ([$]) and and 
hence this must be the pure effect of q- deformation. If we could analyze a more strictly 
g-deformed model, we would observe similar unstable phase transitions in some regions 
of the phase space and in more complicated way. Therefore g-deformed field theory, as 
a formulation of non-archimedian, non-commutative or foam-like structures of spacetime 
fL9f , [20|] , could be related to a certain drastically changing phase structure. Although we 
only studied the D = 2 model, this picture will probably be same even if in higher D 
cases. 

It is still unknown how phase structures depend on details of Q's. To know this, there 
would be no other way to collect many examples. Choosing various functional form of 
Q, it might be interesting to compare how each critical surface differs. For example, for 
the undeformed model, which has another interest of chiral symmetry breaking, Q(s) is 
given by scots in the electric dominant case. In this case, we have to deal with an infinite 



number of singularities of Q. Apart from the singularity problem, various features of 
our oscillatory Q might call us a precaution when studying phase structures of other 
oscillatory charge functions as well as that singular case. 

Finally, we hope that our approximated analyses would be helpful for quantitative 
analyses of strictly g-deformed (quantum group based) models in both contexts. 
Acknowledgments: We thank S. Mukaigawa for comment and information. 
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Figure Captions 

Fig.l: The 2nd order critical line (/i = 0) on £-T plane. S and M mean solitonic (broken) 
and metallic (symmetric) phases respectively. 

Fig. 2: The situation of /i = 0.59. A is a tricritical point. 

Fig. 3: The situation of ji = 0.601. B is a tricritical point. The dashed line means 1st 
order transitions. Two tricritical points Ai and A 2 sit on the boundary between 
1st and 2nd order critical regions. 

Fig. 4: The phase diagram at /i = 0.607. Bi,B 2 , Ai, A 2 are tricritical points. 

Fig. 5: The phase diagram at /i — 0.7. B 2 , A x , A 2 are tricritical points. 

Fig. 6: The phase line from weak ji expansion on plane. C is tricritical. 



